Finite-size effects in microrheology 
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We propose a model to explain finite-size effects in intracellular microrheology observed 
in experiments. The constrained dynamics of the particles in the intracellular medium, 
treated as a viscoelastic medium, is described by means of a diffusion equation in which 
interactions of the particles with the cytoskeleton are modelled by a harmonic force. The 
model reproduces the observed power-law behavior of the mean-square displacement in which 
the exponent depends on the ratio between particle-to-cytoskeleton-network sizes. 
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I. INTRODUCTION 

Transport of particles as vesicles, lipid gran- 
ules or chromosomes through the intracellular 
environment of eukaryotic cells has been subject 
of intense theoretical and experimental research 
in recent yearsi '^i^i^i'^i^i^i'^i^i^i'^^i'^^i^^ In these sys- 
tems, the presence of molecular motors as vi^ell 
as of the cytoskeleton may introduce strong spa- 
tial inhomogeneities and elastic forces which are 
responsible for anomalous transport of the par- 
ticle. 

A manifestation of this anomalous transport, 
observed through video-based and diffusing wave 
spectroscopy methods, is the power law behav- 
ior of the mean-square displacement (MSD) as 
a function of time i^i^i7i?i?i^?i^^i^^i^? These mi- 
crorheology experiments have shown that the 
power-law exponent depends on the ratio be- 
tween the size of the particles to the character- 
istic length of the polymer network constitut- 
ing the cytoskeletont^ i^i^i^ There are two limiting 
cases. The first one in which the radius of the 
particle a is larger than the charateristic length 
of the polymer network ^, o/^ > 1; and the sec- 
ond one in which the radius of the particle is 
smaller than ^, a/^ < 1. 

The results of video-based microrheology ex- 
periments in in vitro F-actin networks shown in 
Fig. 1 are an example of the MSD behavior in 
the second situation. It is important to stress 
that these results are consistent with those ob- 
tained from living cells and in vitro microtubule 



arrays by using also video-based methods, in 
which the MSD of probe particles also shown 
subdiffusion with different exponents)^ 




FIG. 1: MSD vrs time as obtained from the pro- 
posed model for different particle-to-network-size ra- 
tio a/^. Solid lines correspond to the solution of 
Eq.lO with H12|l whereas dotted lines represent the 
power-law given in Eq. 115|l and valid at short 
times. Experimental data (symbols) have been taken 
from Ref. [8] and correspond to particles with ra- 
dius a = 0.25^m and F-actin characteristic network 
sizes: O.lbfim (triangles), O.bbfim (circles), O.SO^m 
(squares) and 0.25^to (diamonds). 

In this article, we give an explanation of 
this effect in the case when the particles have 
a characteristic length smaller than the char- 
acteristic length of the polymer network. We 
propose a diffusion model, taking into account 
the viscoelastic and bounded natures of the 
intracellular medium, to analyze the behavior 



of the MSD through a non-Markovian Brown- 
ian dynamic s '^^i"'^^'^i"'^^i'^'i'i^?i^?i^9i^^i^^ which consid- 
ers particle finite-size effects and confination. 

At short times, om' results lead to a power 
law behavior for the shear modulus G"{uj) as a 
function of the frequency uj in which the expo- 
nent is a function of the particle-to-network-size 
ratio a/^. We also show the relation between the 
measured mean square displacement with an ef- 
fective diffusion coefficient D(t) and an effective 
friction /3(i). It is necessary to stress that these 
coefficients, containing memory effects through 
its dependence on timefiS*ii are related with 
the position relaxation function of the Brown- 
ian particle and, only in the Markovian approx- 
imation, they reduce to the usual diffusion and 
friction coefficients. 

The article is organized as follows. Section 
II will be devoted to introduce the generalized 
diffusion equation for a Brownian particle of fi- 
nite size and the main ingredients of the model. 
In Sec. Ill we apply the model to describe re- 
sults of microrheological experiments in the in- 
tracellular media and in vitro systems. In the 
discussion section we will summarize our main 
results. 



II. THE GENERALIZED DIFFUSION 
MODEL 

For times much larger than the characteris- 
tic time Pq = m/6Trar]s, where a is the radius 
of the particle, m its mass and r]s the viscos- 
ity of the solvent, the dynamics of the Brownian 
particle through a viscoelastic medium can be 
described in terms of a non-Markovian diffusion 
equation of the formpi 



dp 



D{t)V^p-p-'{t)V ■ pF{x 



:i) 



in which the effective diffusion coefficient D{t) is 
related to the effective mobility f3~^{t) through 



D{t) 



ksT 



m 



p-\t), 



(2) 



where ksT is the thermal energy. In Eq. (^, 
p{x, t) is the probability density, and F the force 
accounting for the interactions of the particle 



with the viscoelastic medium. It contains the 
elastic forces exerted by the polymer network 
over the particle and takes into account finite- 
size effects in the dynamics of the particle. 

In the absence of molecular motors the dy- 
namics of the particle is sub diffusive i^*^ In 
this case one may assume, within an effective 
medium approximation, that its motion is cou- 
pled to the medium through a harmonic force 



F{x) 



-ujQex, 



(3) 



where cvq is a characteristic frequency of the 
polymer network. 

The parameter e incorporates the mentioned 
finite-size effects of the particles. In particular, 
it reflects the effects of the stresses exerted by 
the fluid over the surface of the particle during 
its motion though the fluid. Its general form, 

a^/3 • (3, has been calculated in Ref. [20] 



6fcsT 



for the case of a particle moving into a fluid un- 
der flow conditions by taking into account the 
Faxen theore m^^i^^ : the drag force over the par- 



5 

tide is Fdrag (X v — vq with V the velocity of the 

particle, and vq ~ uq + fr^^^o the average of 
the velocity of the fluid vq over the surface of the 
particlei^ For convenience, we will assume that 
e = ToPitQ) with T£) = ma^ (3{tQ) / QkBT the char- 
acteristic diffusion time over the distance a.— 
Therefore, the form of the elastic force © then 
assumes that it originates from the fact that par- 
ticles have finite size and that point particles are 
not affected by the network in the time interval 
we are considering. 

It is important to mention that the non- 
Markovian character of equation (^ is in- 
corporated through the time dependence of 
the effective mobility /3~^(t). This time 
dependence can be shown to arise from a 
description given in terms of a generalized 
Langevin equations j^^i^^i^^i^'^i^^i^^i'^^ or from a 
slow-varying-field approximation of the corre- 
sponding equation obtained from the theory of 
continuous time random walks i^ Also impor- 
tant to mention is the fact that, in the gen- 

eral case, the effective mobility (3 (t) is not the 
usual transport coefficient of the Markovian ap- 
proximation. However, it can be shown that 



it recovers its usual form in the corresponding 
approximation ji£i2ii 

During the motion of the particle through 
the viscoelastic fluid, hydrodynamic interactions 
with the polymer network introduce corrections 
to the mobility of the particle which are simi- 
lar in form to those arising form the presence of 
other particles or from the presence of a wall 
or a liquid-like membranej24i2Si2Si21 For exam- 
ple, in the presence of a wall, these corrections 
depend on the ratio a/h in the form: /3~„;; — 
Pq (1 — 0.625a//i) with h the distance to the 
walli^ It is then plausible to assume that the dis- 
tance of the particle to the "wall" of the "cage" 
formed by the polymer network is of the same 
order as the characteristic length of the polymer 
network ^, and then write in first approximation 



with 



/3-Hi) = /?o 



l-B, 



i 



r'{t), (4) 



where (3 ^(t) is the dimensionless function of 
time accounting for the memory effectsi^ii^iiSiSS 
As already mentioned, the coefficient (l — Bit) 
introduces finite-size effects through the ratio 
a/^ and the parameter Bi takes its value de- 
pending on the nature of the boundaryi ^'^i^^i^^ 



III. CONFINED DIFFUSION OF A 
FINITE SIZED PARTICLE 

Once introduced the general aspects of the 
model, we will analyze the particular case of 
a particle moving in a viscoelastic medium in 
which the effects of constrainng are also present. 
Finally, we will compare our results with exper- 
iments. 

Substituting the harmonic force © , the rela- 
tion for the effective diffusion coefficient (jSJ and 
the effective mobility (jlj), the generalized difu- 
sion equation (^ takes the form 



I -"-('> 



ksT 



m 



VV + e^o V • (xp) 



■ (5) 



A non-stationary solution of Eq. © can be 
obtained by first introducing the dimensionless 
coordinates x = a~^x and the dimensionless 
time i = t^^t, with to = /3q^, and scaling time 



(3-\i')di'. 



(6) 



Assuming a probability distribution of the form: 

Pix,T) = f{T)exp ^-x /A{t) 

solution 

P = 



Po 



\ _ g-Scj^T/aw^ A 2 



one arrives at the 



(7) 



exp 






in which pQ is a normalization constant and we 
have defined the frequency uj"^ = kBT/ma? and, 
according to Eq. @, a is given by 



a 



1-B, 



^ 



(8) 



Eq. (O gives the probability distribution for a 
single particle in the case of constrained diffu- 
sional motionii^ It leads to the following expres- 
sion for the mean square displacement of the par- 
ticle 



(f')(r) 



OiUn-' o 



Pl-^l 



1 + coth 



2^0 



auj' 



T 



(9) 



which can be obtained by using the definition 
{x ) = J X pdx. In Fig. 2, we show a fitt (with 
Eq. Q) of the experimental data of constrained 
diffusion of chromatin inside the nucleus of living 
cells, see caption of the figure for details. 

Experiments measuring the mean square dis- 
placement of Brownian particles in in vitro F- 
actin networks shown a power law behavior char- 
acteristic of subdiffusion in which the exponent 
is a function of the particle-to-network-size ratio 
a/^ (see Fig. 1).-^ To explain the origin of this 
short-time power law behavior of the MSD, we 
will first obtain the explicit dependence of T{t) 
at short times. To achieve this objective, it is 
convenient to analyze the evolution equation for 
the relaxation function x{i) = I x ■ xopdxj^^^ 



d_ 



X{i) 



^0 



r\Mi)- 



(10) 
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FIG. 2: Fitting of experimental data (circles), corre- 
sponding to constrained diffusion in living cells given 
in Fig. 2b of Ref. [12], by means of the MSD ob- 
tained in Eq. Q with lo^/ujo — 4/5, i?2 — 1/6 and 



rameters are given in the text. 



0.02/ir7i^. The values of the pa- 



Eq. ()1U() can be derived by using Eqs. Q and 
Q^i6j20 Then, by taking into account Eqs. © 
and (|1U() one may obtain the relation 



r = In 



X(^'o) 
Xii) 



9 / 9 



(11) 



which can be expressed in a more convenient 
way by expanding x{i) up to first order in t: 
X{i) ~ x(io)(l - Bit) + 0{P), with B2 a param- 
eter characterizing the type of relaxation. The 

resulting expression is r ~ In [l -|- B2i\ ^ " . 
After expanding the logarithm around the unity, 
one finally finds 



Bni°"^T/'^0, 



(12) 



With this expression, the relaxation function is 
found to be 



x(t) ~ x(to)e 



irB2t 



2 / 2 



(13) 



which is a stretched exponential characteristic 
of the relaxation of systems with strong spatial 
inhomogeneities . 

The MSD can now be calculated by using 
the short time approximation of Eq. ((T)), which 
can be obtained by expanding the exponential 



e °"^T up to first order in r, from where one 
obtains 



p{x,t) 



Po 



2lOqUJj^ a "r 



1 xl 



fl2 =.2 



(14) 



Then, from Eqs. H14|) and 1)12(1 one obtains in 
the original variables 



{x^)it)^6DoPo^aB2 
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J. -| auj^/ujQ 



(15) 



where Dq = kBT/mPo is the usual expression 
for the diffusion coefficient in terms of the Stokes 
friction coefficient (Sq. The result clearly shows 
that the finite size of the particles affect the ex- 
pression of the MSD through the factor o/^. 

In Fig.l, we compare our results for the 
MSD with the experimental data reported in 
Ref. [8]. The slopes of the dotted lines were 
obtained by using Eq. ((T3|) and the relation 
a = 1 — Bia/^ with Bi = 2/3 (a value similar to 
that obtained from hydrodynamics^iiSI) Dq ~ 
W^'^cm'^s-'^, l3o ~ 106s-^ LOT ~ 2.25 • lO^s"^ 
and Wo ~ 2 • 10^s~^ for the different values of 
a/S, indicated in the caption. With these val- 
ues of the constants, the experimental magni- 
tude of the MSD were well fitted by assuming 
that a ~ [1 - (2/3)(a/0]exp[-5(a/0]. Finally, 
in figure 1 solid lines represent the fit of the data 
using Eq. Q and thus indicating the effects of 
confinement due to the material surrounding the 
particle. 

At short times, the viscoelastic properties 
of the intracellular medium become manifest 
through the effective diffusion coefficient D(t) 
given in Eq. 0, and which can explicitly be 
obtained by taking the time derivative of Eq. 
((T^ i^?i^7i^? After using the definition of (x^) 
and Eqs. Q, © and (jH)), we have identified 
B2 = a~^ujQ/ujj^ and the explicit expression for 
the dimensionless function P~^{t) which incor- 
porates memory and finite-size effects 



r\t) 



9 I 9 

J. -| auj^/ujQ- 



(16) 



The relation between the effective mobility 
per unit mass (jlj) with the effective viscosity of 



the intracellular medium can be established by 
realizing that P{t) represents a temporal average 
of the for m^^i^^i^° 



m 



C{z)dz, 



(17) 



where C(t) is a memory function related with the 
frequency-dependent effective viscosity through 
the inverse Laplace transform of C('-^) ~ 
?7(a;)ri*22i^ and we have introduced td as a cut- 
off time. After calculating (^{t) from Eqs. (jlGj) 
and (|17j) . the Laplace transform of the result 
leads to the scaling law for the effective viscosity 



r]{uj] 



2 / 2 



(18) 



From this expression, one may calculate the de- 
pendence of the shear modulus G"{lo) = ujri{io), 



as 



G" 



~ UJ 
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(19) 



As in the experiments;^ this expression shows 
that the shear modulus G"{uj) inferred through 
the generalized Stokes-Einstein relation D[lo) 
%Tran(u}) ' depends in general on the ratio a/£,. 
The results here obtained can also be used to 



8.31 



explain the subdiffusive behavior, (x^) 



exp 



t^/\ 



found in experiments with passively diffusing 
particles in living cells »^ The form of the MSD 
with au'^/ul = 3/4, and Eqs. ((111), dl) and 
(|19)) then imply the behaviors [3{t) ~ t^^''^, 
'r]{uj) ~ u>~^''^ and G" ~ u>^'^, in good agreement 
with the results of Refs. [11]. The experiments 
were performed with particles having diameters 
in the range a ~ 2 — Sfim^ Thus, our relation 
a/^ ~ 2/5 implies that the characteristic length 
of the cytoskeleton network is ^ ~ 5 — 7.5^m 
which lies within the same range of values as the 
actin persistence length ^acUn — 5 — 20/umi^ 



IV. DISCUSSION 



We have shown that the observed finite-size 
effects in the anomalous transport of particles in 
the intracellular environment can be explained 
on the grounds of a diffusion model which in- 
corporates the viscoelastic nature of the intra- 
cellular medium, and the characteristic lengths 
associated to the size of the particles and of the 
polymer network constituting the cytoskeleton. 

The proposed diffusion equation incorporates 
memory effects through time dependent coeffi- 
cients which enables one to calculate the mean- 
square displacement from which the role played 
by the mentioned finite-effects effects can be 
analyzed. These effects correspond to the size 
of the particles and the characteristic length of 
the polymer network and of the size of the sys- 
tem, i.e., our description may include the ef- 
fects of confinement. As a consequence, by as- 
suming that the constants take values typical 
of experiments, we have found good agreement 
between our theoretical results and the experi- 
mental data taken from microrheological exper- 
iments with viscoelastic fluids similar to the in- 
tracellular medium. Good agreement has also 
been found in describing confined diffusion of 
chromatin within the nucleus of live cells. 

In summary, the scheme presented can be ap- 
plied to explain the observed subdiffusive behav- 
ior of particles passively diffusing in the cell. 
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